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COMPUTATIONAL MODELING OF VORTEX GENERATORS FOR TURBOMACHINERY 


R. V. Chima 

National Aeronautics and Space Administration 
Glenn Research Center 
Cleveland, Ohio 44135 


ABSTRACT 

In this work computational models were developed and used to 
investigate applications of vortex generators (YGs) to turbomachinery. 
The work was aimed at increasing the efficiency of compressor compo- 
nents designed for the NASA Ultra Efficient Engine Technology 
(UEET) program. Initial calculations were used to investigate the phys- 
ical behavior of VGs. A parametric study of the effects of VG height 
was done using 3-D calculations of isolated VGs. A body force model 
was developed to simulate the effects of VGs without requiring compli- 
cated grids. The model was calibrated using 2-D calculations of the VG 
vanes and was validated using the 3-D results. Then three applications 
of VGs to a compressor rotor and stator were investigated: 1. The 
results of the 3-D calculations were used to simulate the use of small 
casing VGs used to generate rotor preswirl or counterswirl. Computed 
performance maps were used to evaluate the effects of VGs. 2. The 
body force model was used to simulate large part-span splitters on the 
casing ahead of the stator. Computed loss buckets showed the effects of 
the VGs. 3. The body force model was also used to investigate the use 
of tiny VGs on the stator suction surface for controlling secondary 
flows. Near-surface particle traces and exit loss profiles were used to 
evaluate the effects of the VGs. 

NOMENCLATURE 

A c , A vg cell area, vortex generator planform area 
B body force vector 

c chord 

C D lift and drag coefficients 

E, G, H inviscid plus viscous fluxes 

e total energy 

h vortex generator height 

i unit vector 

J Jacobian of transformation from (x, y, z) to (£, r\, q) 

L, D lift, drag 

— ^ 

Fl, Ed lift and drag body forces 

K centrifugal force source term 

P, T, V pressure, temperature, velocity 

q vector of conserved variables 

x, y, z Cartesian coordinates 

x, 0, r cylindrical coordinates 

u , v, w Cartesian velocity components 


u, v 0 , v r 

a 

5 

p 

a 

Q 

co 


cylindrical velocity components 
flow angle 

boundary layer thickness 
computational coordinates 
density 

vortex generator stagger angle 
blade row angular velocity 
loss coefficient 


Subscripts 

0 

1,2 

LE 


stagnation state 
inlet, exit 
leading edge 


Superscripts 

V' relative velocity 

/ unit vector in the direction of the lift 

q q/J - divided by the Jacobian 

INTRODUCTION 

Vortex generators (VGs) are small, vane-like devices that generate 
vortices at their tip. They are commonly used to delay separation on 
aircraft and often used to fix aerodynamic problems that arise during 
testing. In a typical application the vanes are mounted normal to a wing 
or fuselage surface in an array that resembles a cascade of turboma- 
chinery blades. The vanes are staggered at an angle that produces an 
incidence angle with the flow, and the resulting pressure difference 
across the vanes generates vortices at the tip. The height of the VG is 
usually slightly higher than the boundary layer, such that the vortices 
mix high momentum flow from the free stream with low momentum 
flow near the wall. This increase of momentum near the wall allows the 
boundary layer to resist adverse pressure gradients. Although VGs can 
reduce separation, they can produce a measurable increase of drag at 
cruise (Roskam and Lan, 1997.) The reference by Pearcey (1961) gives 
an excellent overview of VG theory and applications. 

Vortex generators may be grouped into two classes: corotating and 
counterrotating. Corotating VG vanes are all staggered in the same 
direction and produce vortices of the same sign. Corotating vortices 
tend to remain at a constant distance from the surface and persist far 
downstream. The vortices interact with each other and with the surface 
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to turn the flow in the direction of the vane stagger. This characteristic 
can be exploited to redirect secondary flows on the surface. Corotating 
vane-type VGs are the most commonly used type of YG and were the 
only kind considered in the present work. 

Counterrotating VG vanes are grouped in pairs with alternating stag- 
ger angles and produce vortices of alternating sign. Depending on the 
blade spacing, the vortex pairs can be made to move towards each other 
and away from the surface, or away from each other and towards the 
surface. This characteristic can be used to control the range of effec- 
tiveness of the vortices but must be carefully assessed for each applica- 
tion. Counterrotating YGs were not considered in this work. 

Many other types of devices have been used as VGs. Vanes, wings, 
ramps, wedges, leading-edge devices, and wall jets have all been used 
as VGs. The reference by Liou et al. (1999) gives an interesting com- 
parison of 12 types of VGs and their effects on heat transfer augmenta- 
tion and skin friction in a square duct. 

In the present work CFD models were used to investigate applica- 
tions of VGs to turbomachinery. The work was aimed at increasing the 
efficiency of compressor components designed for the NASA Ultra 
Efficient Engine Technology (UEET) program. One of the program 
milestones is the design and test of a highly loaded two-stage compres- 
sor (Larosiliere, et al. 2001 and Medd, et al. 2001.) The design point 
pressure ratio is 4.65 at a polytropic efficiency of 92 percent. Two prop- 
erties of VGs were investigated as means for improving the efficiency 
or operating range of the first stage rotor and stator. First, their ability to 
delay separation was investigated as a means for reducing casing 
boundary layer losses in the rotor. And second, their ability to redirect 
secondary flows was investigated as a means for reducing incidence 
losses and spanwise flow migration in the stator. 

Previous use of VGs for Turbomachinery and Internal Flows 

A few researchers have attempted to use VGs in turbomachinery. 
Law, Wennerstrom, and Buzzell (1976) used VGs to energize the cas- 
ing boundary layer ahead of a transonic fan rotor. They measured an 
increase in pressure ratio, a 1-2 point increase in efficiency, and some 
increase in stall margin. They also tried to use small wedge-shaped 
corotating VGs on the suction side of the rotor blades to reduce shock- 
induced separation. In this application they measured a minute increase 
in pressure ratio but lost about one point in efficiency. According to 
Wennerstrom (1987,) “The results achieved with rotor-mounted vortex 
generators were less than rewarding.” 

McCormick (1993) also used low-profile VGs to reduce shock- 
induced separation on a turbulent flat plate boundary layer. While the 
extent of the separation was reduced, he found that the mass-averaged 
total pressure loss increased. Gammerdinger (1995) found similar 
results in a transonic fan-blade cascade. 

At high altitude flight conditions the Reynolds numbers on low pres- 
sure (LP) turbine blades is very low and separation is a common prob- 
lem. Bons, Sondergaard, and Rivir (2000) used both steady and pulsed 
VG jets to control separation on LP turbine blades in a linear cascade. 
They completely eliminated separation and reduced wake profile loss 
by 50 percent. 

DeCook, King, and Elrod (1993) used crenulated (notched) trailing 
edges on a compressor cascade to increase wake mixing. They mea- 


sured a decrease in total pressure loss of 10-20 percent and a slight 
decrease in flow turning. Kuethe (1972) attached unusual half-cylinder 
shaped VGs near the trailing edges of a flat plate and an airfoil to 
improve wake properties associated with noise. He was able to reduce 
the wake velocity deficit and suppress vortex shedding. Similar results 
were found for a transonic fan blade cascade tested by the acoustics 
branch at NASA Glenn (Lewis) Research Center around 1994 (unpub- 
lished.) However, preliminary tests of Kuethe- type VGs affixed to a 
rotating fan showed little effect on noise generation. 

VGs have been used to control secondary flows in inlet ducts. In a 
collaborative effort between NASA Glenn Research Center and the 
Defence Research Agency, Anderson and Gibb (1998) used VGs to 
control vortex liftoff (separation) in an S-shaped inlet duct. Experi- 
ments showed large reductions in engine face distortion and unsteadi- 
ness. In subsequent work Anderson, et al. (1999) used Design of 
Experiments techniques with a computational fluid dynamics (CFD) 
model to optimize VG geometry for an inlet S-duct. Their results 
showed that very small “MEMS scale” devices, could be used to mini- 
mize distortion. Their devices were reported to be on the order of the 
boundary layer momentum thickness in height, although very long 
devices appeared to be the most effective. Also at NASA GRC, Wendt 
and Dudek (1998) combined CFD modeling and experiments to 
arrange VGs in a subsonic bifurcated diffuser. They were able to reduce 
engine face distortion by 59 percent with only a 0.3 percent reduction 
in total pressure recovery. In all of this work CFD modeling was 
heavily used to optimize VG geometry and placement (Bender, et al. 
1999.) The successful use of VGs for controlling secondary flows 
in inlet ducts motivated many of the present applications to 
turbomachin ery. 

Objectives of the Present Work 

In the present work isolated VGs were modeled computationally 
using direct simulation with the Swift code (Chima, 1996.) A paramet- 
ric study was performed to determine the effects of VG height on total 
pressure mixing and flow turning. Next a body force model was devel- 
oped for Swift that allowed the effects of VGs to be simulated without 
detailed gridding of the VGs. Then three applications of VGs to turbo- 
machinery were studied: 

1 Use of VGs for rotor preswirl. The direct simulation results were 
used to provide mixed-out inlet conditions to the UEET first stage 
rotor to simulate the effects of VGs on the casing. Both prewsirl and 
counterswirl orientations were investigated. Computed performance 
maps were used to evaluate the effects of the VGs on pressure ratio, 
efficiency, and stall range. 

2 Use of VGs as part-span splitters. The body force model was used 
to simulate part- span splitters on the casing ahead the UEET first 
stage stator. Computed loss buckets were used to evaluate the effects 
of the VGs. 

3 Use of blade surface VGs for secondary flow control. The body 
force model was also used to simulate a spanwise array of VGs 
along the suction surface of the stator. Many configurations were 
examined to try to reduce secondary flows. Secondary flow patterns 
and spanwise loss profiles were used to evaluate the effects of the 
VGs. 
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Figure 1. — Computed particle traces and PO contours, 

VG height = 1.3 5. 

CFD CODES 
RVCQ3D 

The quasi-three-dimensional turbomachinery analysis code 
RVCQ3D developed by Chima (1987, 1995) was used to generate lift 
and drag coefficients for the VG model. The code solves the thin-layer 
Navier-Stokes equations in finite-difference form. Radius change and 
variable stream sheet thickness can be modeled but were not used here. 
Turbulence effects were modeled using the Baldwin-Lomax model, 
including the original transition model. The flow equations were solved 
using an explicit Runge-Kutta scheme. A spatially varying time step 
and implicit residual smoothing were used to accelerate convergence. 
Solutions were run on a C-shaped grid with 191 x 45 points. Typical 
solution times were about one minute per case on an SGI Octane work- 
station. 

Swift 

All three-dimensional solutions were run using the Swift multiblock 
turbomachinery analysis code developed by Chima (1996.) Swift 
solves the three-dimensional Navier-Stokes equations using algorithms 
similar to RVCQ3D. Multiblock capabilities allow the tip clearance 
region over a blade to be computed directly, so isolated VGs were mod- 
eled as stators with a very large clearance. The Baldwin-Lomax turbu- 
lence model was used for all rotor cases based on considerable past 
experience. The k-co model was used for all stator/VG applications and 
isolated VG cases, where resolution of discrete vortices was expected 
to challenge the algebraic model. Most 3-D computations were run on 
the Cray C90 at NASA Ames Research Center (von neumann) and took 
about two hours per case. Swift autotasks well so that most cases took 
less than 30 minutes wall clock time. 

DIRECT CALCULATION OF ISOLATED VGS 

Isolated VG vanes were analyzed using Swift to examine the effect 
that the VGs had on the flowfield and to investigate the effects of VG 
height on their performance. Although VGs are often made from flat 
plates, it was felt that more sophisticated vane shapes might be used in 
a turbomachinery application. Hence, a NACA 64-006 profile with a 
rounded trailing edge was used. The vanes were sized to fit on the cas- 
ing in front of the UEET compressor rotor in an arrangement like that 
used by Law, et al. (1976.) Using five vanes per rotor passage plus one 



Figure 2. — Computed particle traces and PO contours, 

VG height = 5/4. 

to reduce resonances gave 131 vanes. Common practice for VGs is to 
use a solidity of one and a stagger angle of 20 degrees. This gave an 
aerodynamic chord of 12.7 mm (0.5 in.) 

Computational grids were generated for the vanes. The grids had two 
blocks. A C-shaped block was used around the vanes with 191 points in 
the C-direction and 45 points from vane to mid-passage. Sixty three 
points were used span wise, with 15-31 points along the vane depend- 
ing on height. An O- shaped block filled the hole above the tip of the 
vanes. The total grid size was about 554K points. The grid extended 4.4 
VG heights above the VGs and 2 chords downstream. 

Boundary conditions for Swift were specified as follows: At the inlet 
boundary total temperature and flow angles were specified as constant, 
and a total pressure profile was specified to match an equilibrium turbu- 
lent wall-wake velocity profile with a boundary layer thickness 
5 = 1.9 mm (0.075 in.) At the exit boundary the static pressure was 
set to give a free stream Mach number M = 0.53 . 

Four vane heights were analyzed, h - 1.35, 5/2, 5/4, and 5/8 , 
where 5 = 1.9mm (0.1 in.) . The largest vane height, h = 1.3 was 
five times the rotor tip clearance. The smallest vane height, h - 5/8 
was on the order of the boundary layer displacement thickness. Since 
the vanes were immersed in a turbulent boundary layer, the flow over 
them was modeled as fully turbulent using the k-co model. 

Figure 1 shows the computed solution for h = 1.35 . Particle traces 
colored by total pressure leaving the tip of the vane show the structure 
and path of the tip vortices. The core of the vortices has low total pres- 
sure but the outer portion of the vortices clearly transports high total 
pressure flow from above the vanes to the wall. Total pressure contours 
on the exit plane show how boundary layer height is reduced on the 
pressure side of the passage and increased on the suction side. Figure 2 
shows similar results for h = 5/4 . Here the oncoming boundary layer 
flow has much lower momentum and the vortex strength is greatly 
reduced. 

The flow at the exit planes was energy-averaged circumferentially 
and the resulting spanwise profiles of total pressure and flow angle are 
shown in figs. 3 and 4 respectively. The unmodified boundary layer 
profile is shown for reference. The total pressure profiles in fig. 3 show 
that the two largest VGs transport core flow to the wall. The profiles 
near the wall are similar for h - 1.35 and h - 5/2 , but the VGs with 
h = 5/2 have less impact on the core flow. The two smallest VGs 
have little effect on the total pressure profiles. 
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Figure 3.— Spanwise P0 profiles two chords downstream Figure 6.— Body force nomenclature for VG model, 
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Figure 4. — Spanwise turning angle profiles two chords 
downstream of VGs. 
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Figure 5. — Effect of VG height on P0 loss. 


The turning angle profiles in fig. 4 show that the turning decreases 
quickly with vane height. Overall total pressure losses calculated 
between the grid inlet and exit are shown in fig. 5. The two smallest 
VGs show little additional loss over the unmodified boundary layer. 
VGs with h = 5/2 increase the baseline loss about 60 percent, but 


VGs with h = 1.35 nearly quadruple the loss. The results of this 
parametric study suggest that h = 5/2 is about optimal for mixing the 
flow with minimum loss, but that tall VGs must be used if the intent is 
to turn the flow. 

VG MODEL 

Formulation 

A computational model of vortex generators was developed and 
implemented in the Swift code. The model uses body force terms to 
produce the effects of VGs on the flow without the difficulties of grid- 
ding and solving the flow around the VGs directly. The model was 
based on the one developed by Bender, et al. (1999.) Bender’s model 
used a lift force normal to the flow to turn the flow. The lift force was 
proportional to the difference between the desired flow angle and the 
actual flow angle. The constant of proportionality was shown to have an 
“asymptotic mode” in which the strength of the induced vortex was 
independent of the value of the constant. No drag force was included so 
that the model produced no losses. 

The present model includes a lift body force that turns the flow and a 
drag body force that produces losses. The body forces were assumed to 
be proportional to standard lift and drag coefficients that are functions 
of the incidence angle. The lift and drag coefficients were calibrated 
using two-dimensional RVCQ3D calculations of a cascade of VGs. The 
model described below applies to VGs on the hub or casing of a turbo- 
machine. Similar equations may be derived for VGs on the blade sur- 
faces with a simple change of variables. 

Figure 6 shows a cascade of VGs stacked along a radial line and set 
at a stagger angle o with respect to the x-axis. The flow direction is at 
an incidence angle a with respect to the ^tagger angle. The aerody- 
namic forces on the blades are a lift force L normal to the flow direc- 
tion and a dra^ force D in the flow direction. The body forces on the 
flow F L and F D are equal and opposite to L and D . 

For the Swift code the Navier-Stokes equations are written in Carte- 
sian coordinate system rotating around the x-axis with angular velocity 
Q . The equations are transformed to a general body-fitted coordinate 
system using standard finite- difference techniques (Chima, 1992) and 
may be written as: 
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d t q + d e E + O^G + d q H = K + B (1) 

where q = (p, p u, pv, pw, e)/J is the vector of conserved vari- 
ables, E, G, and H are vectors containing the inviscid and viscous 
fluxes, and k contains centrifugal force terms for the y- and z-momen- 
tum equations. Here the caret notation denotes division by the Jacobian 
of the coordinate system transformation J . The body force term B is 
given by 


0 

(Fl + K)x 

(p L +K) y 

ra(F L + F* D ) 0 


( 2 ) 


where F 0 = ( E y z-F z y)/r and 


~f l = ijplvl 2 C L A c ^-l 


C l A c 


K = l -p\vYc D A c ^-d 


(3) 


reduce discontinuities in the finite-difference solution. The model may 
also be applied over n consecutive surfaces to generate the overall 
turning and loss produced by the VGs without resolving discrete vorti- 
ces. In this case - 1 /n. 

The lift force is assumed to be normal to both the local relative flow 
direction and the spanwise (radial) direction. In this way it can turn the 
flow without producing losses. Using carets to denote unit vectors, the 
unit vector in the relative flow direction is 
v = (ui x + Vq'/q + v r i r )/\ V'\ . The unit spanwise vector is 
s - 0 i x + O/ 0 + s r i r , where s r = 1 for YGs on the hub or s r - -1 
for VGs on the casing. The unit vector in the direction of the lift is 
given by: 


l = V'xs 
= |Pj( v e l x~ ul e) 


( 6 ) 


where v 0 ' = (v'z-w'y)/r and i 0 = (i y z — i z y)/r. 

The drag force is assumed to be parallel but opposed to the relative 
flow direction so that it can produce losses without turning the flow. 
The unit vector in the direction of the drag is given by: 


d =-(ui x + Vi + w'i z ) /\V'\ 


(7) 


Calibration 


Here 


ip| V'| 2 = ^p[w 2 + v' 2 + w' 2 ] (4) 

is the dynamic head in the relative frame of reference and 
v' = v-Qy, w' = w + £lz are relative velocity components. 
C L and C D are lift and drag coefficients, A c and A vg are areas, and 
/ and d are unit direction vectors described below. 

The energy equation includes a term for the work done on the fluid 
by rotating VGs, for example_by VGs mounted on a rotating hub. It 
could also include a term like V' • F D , but the model assumes that this 
term exactly cancels the entropy rise term due to the drag forces. In this 
way the drag model does not change the total enthalpy of the flow (see 
Marble, 1964.) 

The VG model is applied on grid surfaces that approximate the VG 
planform. For VGs on the hub or casing, the model is applied on a 
(8, q) surface, and the r \ -coordinate is approximately normal to that 
surface. Then A c is the cell area in the direction of the planform, and is 
given by: 


= ^9 

c rJ cos a 


(5) 


The lift and drag coefficients C L and C D were calculated for a 2-D 
cascade of VG blades using the RVCQ3D code (Chima, 1996.) The 
blade shapes and flow conditions were the same as those used for the 3- 
D simulations described earlier, i.e., solidity = 1, stagger angle = 20 
degrees, and M = 0.53 . Several incidence angles were run. The pres- 
sure forces and shear forces in the x- and y-directions were integrated 
around the blades and output, along with the inlet dynamic head and 
inlet and exit flow angles. 

An average flow angle a was calculated and the forces were 
resolved into components tangent and normal to the average flow direc- 
tion, similar to the procedure given by Dixon, (1975.) 

a = atanj^tanocj 4 tana 2 )J 

L = -/^sina + /^cosa 
D = /^cosa + f y sina 

C = L (8) 

\pW'\ 2 ch 


\p\V'l 2 ch 


where T| 0 = (r|^ e + r| z z e ) = (-r^z + ri^) . 

A l , f ,/y / 4 ( , is the ratio of the planform area of the real VGs to the 
total cell area on which the model is applied. The model is usually 
applied on grid surfaces with the same total area as the real VG plan- 
form. To generate discrete vortices the model is applied over two neigh- 
boring surfaces, and A c = 1/2 . Two surfaces are used to 


where c is the chord and the height h - 1 . The results are plotted in 
fig. 7 as a function of incidence angle a, which can be approximated as: 
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Figure 7. — Calculated lift and drag curves for VG cascade, 
a =20 deg. 



Figure 8. — Particle traces and P0 contours computed 
using VG model, VG height = 1.3 5. 


800 



Figure 9. — Comparison of modeled and resolved P0 
profiles. 



Figure 10. — Comparison of modeled and resolved turning 
angle profiles. 


a ~ sin a = V' LE ■ h 

ft = -sina/^H- cosa/ 0 + 0z' r (9) 

a ~ (-u LE sina + v' LE cosa) /\V'\ 

C L is linear for -10° < a <40°, and C D is linear with C E . 
C L and C D were curve fit using: 

C T ~ 2.356a 

2 («>) 

C D = 0.019 + 0.1637C£ 

These coefficients are only valid for VGs with a solidity of 1 and a 
stagger angle of 20 degrees in turbulent flow. They should be nearly 
independent of Mach number as long as the flow remains subsonic. 


Validation 

The VG model was used to recompute the solution for the isolated 
VGs with h - 1.35 shown in fig. 1. For these calculations the 
191 x 45 point C-grid was replaced with a uniform 97 x 27 point H- 
grid, but spanwise clustering was unchanged to maintain the boundary 
layer resolution. The body force model was applied on two grid planes 
near the center of the grid. 

Particle traces and exit total pressure contours computed with the VG 
model are shown in fig. 8 and may be compared with the direct calcula- 
tions shown in fig. 1, which has the same contour levels. The modeled 
vortex lacks the detailed structure of the direct solution and does not 
pick up the low pressure core of the computed vortex, but the effects on 
the exit total pressure contours are quite similar. Figure 9 compares the 
spanwise total pressure profiles from the model with that from the 
direct calculations. The model overpredicts the mixing but does match 
the wall shear and overall losses closely. Figure 10 shows that the 
model matches the turning almost exactly. 
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Figure 11. — Casing VGs ahead of UEET rotor 1 tip section. 



Figure 12. — Effect of pre and counterswirl VGs on UEET 
rotor performance. 

APPLICATIONS TO TURBOMACHINERY 
Use of VGs for rotor preswirl 

Law, et al. (1976) used VGs to energize the casing boundary layer 
ahead of a transonic fan rotor, reasoning that a healthier boundary layer 
would delay stall and possibly increase efficiency. They used corotating 
VGs staggered in the same direction as the rotor, based on an argument 
that gyroscopic forces induced by blade-to-blade pressure forces would 
keep the vortices close to the casing. They sized their VGs based on 
intuition and on recommendations in Pearcey (1961,) using a height of 
2.4 times the rotor clearance, a stagger angle of 20 degrees, and about 
4.8 VGs per rotor passage. They measured an increase in pressure ratio, 
a 1-2 point increase in efficiency, and some increase in stall margin. 

A similar application was considered here by modeling an array of 
VGs on the casing ahead of the UEET rotor 1, shown at the tip in fig. 
11. The VGs were modeled by using the results of the direct calcula- 
tions of isolated VGs as mixed-out inlet conditions to the UEET first 


stage rotor. The VG height was taken as 
h = 1.35 = 2.54mm (0.1 in) . This height was five times the tip clear- 
ance or about two percent of the rotor span. 

For the rotor calculations, spanwise profiles of P 0 , T 0 , and v 0 were 
specified at the inlet boundary. Over most of the span these quantities 
were obtained from a previous analysis of an inlet guide vane. Over the 
last three percent of the span these quantities were obtained from the 
isolated VG calculations. By changing the sign on v 0 both counter- 
swirl and preswirl VGs could be simulated. Counterswirl VGs were 
staggered like the rotor blades (but much less) and turned the flow 
against the rotor rotation. Preswirl VGs were staggered opposite to the 
rotor blades and turned the flow in the direction of rotor rotation. 

The inlet axial and radial velocity components of the profile were 
determined as part of the solution. The casing static pressure was pre- 
scribed at the rotor exit and radial equilibrium was used to determine 
the nominal spanwise pressure distribution. A characteristic boundary 
condition allowed the exit static pressure to vary tangentially about the 
nominal values. 

A C-type grid was used around the rotor with 209 points around the 
blade, 46 points from blade to mid-passage, and 61 points spanwise. An 
O-type grid was used in the tip clearance region, with 9 points between 
the blade tip and the casing. The total grid size was 607,778 points. 
Operating maps were calculated by changing the exit casing static pres- 
sure. Solutions near maximum flow converged to constant mass flow 
and pressure ratio in 2-3000 iterations. Stall points were determined by 
restarting from a solution at a higher flow and increasing the exit pres- 
sure slightly until the solutions either diverged or the mass flow fell 
well below the expected stall point. The stall point for the case with 
counterswirl VGs was run 19,000 iterations to assure convergence. The 
numerical stall points are known to within 0.23 kg/sec (0.5 lbm/sec,) 
but it is not known how well the numerical stall points corresponds to 
the physical stall points. 

Calculated performance maps of pressure ratio, temperature ratio, 
and adiabatic efficiency are shown in fig. 12. The solid line represents 
the baseline rotor without VGs. Counter swirl VGs (O’s) turn the inlet 
flow in a direction opposite to the blade rotation, thereby loading up the 
tip sections and increasing the pressure and temperature ratios accord- 
ingly. This is consistent with the experiments with counterswirl VGs 
described by Law, et al. (1976.) However, the experiments showed an 
increase in mass flow range, while the calculations show a 12 percent 
decrease in flow range. Preswirl VGs (A ’s) unloaded the tip and 
decreased the pressure and temperature ratios, but increased the flow 
range by nearly 15 percent. Both of the computed VG configurations 
showed a negligible decrease in adiabatic efficiency. 

The present results suggest that small preswirl VGs could have a 
large impact on stall range. However; there are several limitations to 
these calculations. First, prediction of stall points is notoriously unreli- 
able. Second, the present results are steady, but the actual rotor would 
see unsteady loads due to the VGs that could be detrimental to perfor- 
mance. Finally, in the present calculations the beneficial effects of the 
VGs on the casing boundary layer are applied at the rotor inlet bound- 
ary. In reality discrete vortices would convect into the moving rotor and 
could influence the casing boundary layer well downstream. Thus it is 
clear that either experiments or detailed unsteady calculations would be 
required to confirm these results. 
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Figure 13. — Left: tip section of 3-D UEET stator 1. 
Right: VGs + flank-milled stator 



Figure 14. — Particle traces from part-span splitter VGs. 


Use of VGs as part-span splitters 

In a second application, VGs were used as part span splitters to con- 
trol incidence effects on the casing ahead of a redesigned UEET stator. 
These stators are highly three-dimensional, with bow and sweep. The 
tip sections (fig. 13 left) are highly cambered to deal with the over- 
turned flow leaving the rotor. In this application it was proposed that 
the 3-D stator blade could be replaced with a straight (or flank-milled) 
blade (fig. 13 right) that would be easier to build. The straight blade 
would also have less surface area than the 3-D blade, and hopefully 
lower skin-friction losses as well. An array of VGs was to be used on 
the casing to guide the overturned flow from the rotor into the stator 
passages (fig. 13 right.) 

A straight stator was designed with a constant camber line and a lin- 
ear thickness variation along its span. It had the same chord as the 3-D 
blade at the hub and the same turning as the 3-D blade at mid span. 
Three VGs were equally spaced across each stator passage, giving a 
spacing of 10.2 mm (0.4 in.) The solidity was set to one. The VGs were 
staggered at the stator inlet angle of 43 degrees, giving an axial chord 
of 7.5 mm (0.29 in.) The VG height was set to 10 percent of the stator 
span (7.6 mm or 0.3 in,) corresponding to the region of highest rotor 
overturning. 

A computational grid was generated for the stator with 
173 x 46 x 51 = 405, 858 points. The stators were cantilevered and 
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Figure 15. — Effect of part-span splitter VGs on UEET 
stator losses. 


the hub clearance effects were modeled using a simple periodic bound- 
ary condition applied over seven points. This clearance model is exact 
for infinitely thin blades and works quite well for most compressor 
geometries (see Chima, 1996 for details.) 

The VGs were computed using the body force model applied on a 
separate H-grid just ahead of the stator. The lift and drag coefficients 
given in equation (10) were calculated for VGs staggered at 20 degrees. 
For the present stagger angle of 43 degrees the lift and drag coefficients 
were recalculated as: 


C T ~ 3.546a 

, ( 11 ) 

C D » 0.0148 + 0.0607 

Figure 14 shows the arrangement of the VGs ahead of the stator 
blades. Computed particle traces show the tip vortices. In this applica- 
tion the VGs were mainly used to turn the flow, and the vortices were 
not expected to affect the casing boundary layer appreciably. 

The rotor exit P Q , T 0 , and v 0 profiles were specified as inlet profiles 
to the stator. To simulate off-design operation, the rotor exit velocity 
was rotated tangentially ±6 degrees and new v e profiles were gener- 
ated. These profiles were used to generate the loss buckets shown in fig. 
15. The baseline 3-D stator (O’s) has low losses for incidence angles 
between -2 and +6 degrees. Losses increase considerably for angles 
less than -2 degrees. The straight blade alone (D's ) has only slightly 
higher losses than the 3-D blade. The straight blade with VGs ( V's) is 
slightly worse that the straight blade alone, except at the highest inci- 
dence. 

Examination of the solutions for the 3-D baseline blades revealed 
that the high losses at negative incidence were due to high Mach num- 
bers near the hub, and that the blades were fairly insensitive to inci- 
dence effects at the tip. Since the straight blade matched the 3-D blade 
at the hub, its performance was similar to the 3-D blade. In this applica- 
tion the part- span splitter VGs had little effect except to add their own 
losses to the flow. In other applications part-span splitter VGs might 
still prove to be a viable concept. 
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Figure 16. — Pressure contours and secondary flows on 
the suction side of UEET stator 1. 



Figure 17. — Reduction in secondary flows using suction 
surface VGs. 


Use OF blade surface VGs for secondary flow control 

Anderson, et al. (1998, 1999) used VGs to control secondary flows in 
S- shaped inlet ducts. A similar application was considered here by 
using VGs to control secondary flows on the suction side of the UEET 
stator. CFD solutions of these stators show strong secondary flows on 
the suction side. Figure 16 shows suction side pressure contours and 
near-surface flow patterns for a case at design incidence. Secondary 
flows migrate up from the hub and down from the casing, merging near 
mid span. At other incidence angles the secondary flows merge into 
small separation bubbles. In this application spanwise arrays of VGs 
were used to redirect the secondary flows. It was hoped that this would 
reduce losses by shortening the streamlines and reducing the size of 
separation bubbles. 

VGs were arranged spanwise on the suction side of the stators in var- 
ious configurations like the one shown in fig. 16. The VGs had a chord 
of about 6.35 mm (0.25 in) and a solidity of one to give 8-10 VGs 
along the span. The VGs were staggered +20 degrees. Suction surface 
boundary layer heights were estimated from earlier solutions, and VG 
heights were chosen with h = 8/2 = 0.71 mm (0.028 in.) The body 
force model was applied over a region of the grid to simulate the over- 
all effects of the VGs without resolving discrete vortices. The grid size 
was the same as in the previous section. 



Figure 18. — Effect of suction surface VGs on stator 
spanwise los profile. 

Several configurations of VGs were considered. The best results are 
shown in fig. 17. Here eight VGs were arranged spanwise between 
approximately 60 and 70 percent chord. The upper four were staggered 
upwards and the lower four were unstaggered. The resulting surface 
pressures were modified slightly (compare figs. 16 and 17) and the sec- 
ondary flows were successfully reduced, as shown by the nearly 
straight surface streamlines. 

The solutions were energy-averaged circumferentially and the span- 
wise profiles of loss coefficient co = (P Q1 - P 0 2 )/(P 0 i “ ^l) are com_ 
pared in fig. 18. The apparently negative loss near the tip results from 
spanwise redistribution of the flow. Inside the stator passage high loss 
flow from the tip of rotor 1 exchanges places with better flow from 
lower along the span, such that the losses appear negative near the tip. 
With VGs the core loss decreased about a point over much of the mid- 
span region, but the loss increased near the endwalls. The net loss 
increased by 0.2 points with VGs. 

At -6 degrees of incidence the secondary flows were stronger and 
harder to control with VGs. Many VG configurations were examined 
but the results always showed the same trend: it was possible to reduce 
the secondary flows with VGs but the loss always showed up some- 
where else, and the net loss was always slightly worse with VGs. 

CONCLUSIONS 

In this work CFD models were used to investigate applications of 
vortex generators (VGs) to turbomachinery. The work was aimed at 
increasing the efficiency of compressor components designed for the 
NASA Ultra Efficient Engine Technology (UEET) program. Initial cal- 
culations were used to understand the physical behavior of VGs and to 
develop a simple CFD model of VGs. Then three applications of VGs 
to the first stage of the UEET compressor were investigated. The fol- 
lowing conclusions were made: 

1 A parametric CFD study of corotating VGs suggested that the opti- 
mal VG height for mixing the flow with minimum loss is about half 
the boundary layer thickness, but that taller VGs should be used to 
turn the flow. 

2 A body force model was developed that allows the effects of VGs to 
be added to turbomachinery calculations without gridding the VGs 
directly. The model was calibrated using lift and drag coefficients 
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calculated with a 2-D CFD code (RVCQ3D,) and validated against 
3-D calculations of fully gridded VGs. The VG model reproduced 
the wall shear, overall total pressure loss, and turning of the fully- 
gridded calculations correctly, although mixing was overpredicted. 

3 The results of the parametric study were used to simulate VGs on 
the casing ahead of the UEET first stage compressor rotor. The com- 
puted VG exit profiles were averaged circumferentially and used as 
mixed-out inlet conditions near the tip of the rotor. Performance 
characteristics were calculated for the baseline rotor alone and with 
counter swirl and preswirl VGs. The counterswirl orientation gave a 
slight increase in total pressure rise but reduced the stall margin. The 
preswirl orientation reduced the pressure rise slightly, but increased 
the stall margin. Neither orientation changed the efficiency signifi- 
cantly. The results suggest that small preswirl VGs (about two per- 
cent span high) could be used to increase the stall margin of 
compressors with little impact on performance. However, detailed 
unsteady calculations or experiments would be needed to confirm 
these results. 

4 The body force model was used to simulate part- span splitters on the 
casing ahead of a stator. Here the bowed and swept UEET first stage 
stator was replaced with a straight blade plus an array of tall VGs 
(10 percent span high) on the casing to deal with overturned flow 
from the upstream rotor. Loss buckets were calculated for the base- 
line blade and the straight blade with and without VGs. The results 
showed that this stator was relatively insensitive to incidence at the 
tip, so that VGs did not improve the performance of the straight 
blade. Part- span splitter VGs might still prove to be useful in other 
applications. 

5 The body force model was also used to control secondary flows on 
the suction side of the UEET stator. Here an array of very small VGs 
(less than 1 mm high) were used to turn spanwise secondary flows 
back towards the endwalls. The results showed that secondary flows 
could be nearly eliminated at low incidence and reduced at high 
incidence. However, reducing secondary flows did not decrease the 
overall loss in the stator. The concept might still be useful for elimi- 
nating massive separations caused by small secondary flows. 
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